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1. Introduction 

In this paper we construct a finite range decomposition for a family of translation invariant 
gradient Gaussian fields on Z*^ (c? > 2) which depends real-analytically on the quadratic from 
that defines the Gaussian field. More precisely, we consider a large torus (z / L^zY and obtain 
a finite range decomposition with estimates that do not depend on N . Equivalently, we show 
that the discrete Greens function Ca of the (elliptic) translation invariant difference operator 
A = V*AV can be written as a sum Ca = CA,k of positive kernels CA,k which are supported in 
cubes of size ~ with natural estimates for their discrete derivatives V°'CA,k (see Theorem 12. ip 
as well as for their derivatives with respect to A (see Theorem 12. 2p . 

To put this into perspective recall that an R™- valued Gaussian field ,^ on Z*^ (with vanishing 
expectation, e(^(x)) = 0) is said to have range M if the correlation matrices E [,^''(x).^*(y)] , 
r,s = 1, . . . ,m, vanish whenever \x — y\ > M. In the following we consider only translation 
invariant Gaussian fields. We say that Gaussian fields form a finite range decomposition 
of ^ if ^ = X^fcCfc and C,k has range ~ L'^, where L > 2 is an integer. The existence of such 
a decomposition is equivalent to a decomposition of the correlation matrices C^(x, y)''''' := 
E [^^(x)^^(?/)] as a sum of positive (semi-) definite matrix valued kernels CA,k with range ~ L'^, 
i-e., E.,y Er,s CXfc(x, y)Cm'{x) > and C^,fc(x, y) = ii \x - y\ > LK 

We are interested in gradient Gaussian fields, i.e., Gaussian fields with a-algebra deter- 
mined by the gradients V.^. Such fields arise naturally e.g. in problems in elasticity where 
only the difference of values is relevant for the energy. In this case we seek a decomposition 
into gradient Gaussian fields such that the gradient-gradient correlation E[Vi^^{x)VjC,^{y)~\ = 
ViV*E[^^{x)^^{y)~\ vanishes for \x — y\ > L^. Gradient Gaussian fields are more subtle to han- 
dle since they exhibit long-range correlations (the gradient-gradient correlation of the original 
field typically has only algebraic decay \x — y\~'^ with the critical exponent —d). In the language 
of quantum field theory gradient Gaussian fields are thus often referred to as massless fields. 

Decomposition into a sum of positive definite operators has been discussed in [8] where a 
radial function is written as a weighted integral of tent functions. In [5], finite range decompo- 
sitions of the resolvent of the Laplacian (a — A)~^, with a > 0, have been obtained both for the 
usual Laplacian and for finite difference Laplacian on the simple cubic lattice Z*^. In [1] these 
results are extended and generalised by providing sufficient conditions for a positive definite 
function to admit decomposition into a sum of positive functions that are compactly supported 
within disks of increasing diameters More precisely, the authors of [1] consider positive 

definite bilinear forms on and prove that finite range decompositions do exist when the 
bilinear form is dual to a bilinear form ip ^ j \Bif{x)\'^ dx where 5 is a vector valued partial 
differential operator satisfying some regularity conditions. 

The main novelty of our paper is twofold. First we extend the finite range decomposition 
for the discrete Laplacian to a situation where no maximum principle is available (even in the 
scalar case there is no discrete maximum principle for general elliptic difference operator V*AV 
with constant coefficients). This can be seen as an adaptation of [1] to the discrete setting. 
Secondly, we show that the finite range decomposition can be chosen so that the kernels CA,k 
depend analytically on A as long as A is positive definite. 

Our main motivation is the renormalization group (RG) approach to problems in statistical 
mechanics, following the longstanding research programme of Brydges and Yau [6], and in 
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particular the recent work of Brydges [3], both inspired by the work of K.G. Wilson [TT]. The 
goal is to get good control of the expectations e{K) of nonlinear functions that depends on a 
gradient Gaussian field in a large region L C Z'^ of the integer lattice. The size of L and the 
long range correlations in C, make it difficult to obtain accurate estimates on the expectation 
e{K). In [2] we show that such control can nonetheless be obtained in many interesting cases 
using the RG approach. One key difference with the earlier work of Brydges and others is the 
necessity to drop the assumption of isotropy. Hence the relevant quadratic term is a general 
finite difference operator V*v4V and reduces no longer to a multiple of the discrete Laplacian. 
For this reason we need a finite-range decomposition for general (elliptic) operators A and, in 
addition, we need to control derivatives of the finite range decomposition with respect to A. 

In Section [2] we introduce the setting of gradient fields and the relevant Greens functions. 
Our two main results are given in Theorem 12.11 and Theorem 12.21 The existence of the finite 
range decomposition is proved in Section [3] where we adapt and extend the methods in [1] to 
our setting. The regularity estimate for a fixed A is established in Section |H Real-analytic 
dependence on A is proved in Section [5l 

Acknowledgments. The research of S.A. was supported by EPSRC grant EP/I003746/1, 
R.K. by the grants GACR 201-09-1931, 201/12/2613, and MSM 0021620845, and S.M. by the 
DFG Research group 718 "Analysis and Stochastics in Complex Physical Systems" . 

2. Notation and main results 

We are interested in gradient Gaussian fields on bounded domains in Z*^. For that let L > 3 
be a fixed odd integer and consider for any integer N the space 

Vn = W : Z'^ ^R"^; ip{x + z) = ip{x) for all z G {L^Z^} = (r™)'"' 
of functions on the torus Tjv := {z/L^zY equipped with with the scalar product 

= ^{V{x),i){x))^m. (2.1) 

Notice that a function on T^v can be identified with an L^-periodic function on Z'^. We will 
later denote the corresponding space of C "^-valued function, equipped with the usual scalar 
product in the same way. 

We consider two distances on Z"^: p{x,y) = mi{\x — y + z\: z G {L^zY} and Poo{x,y) = 
inf{|x — y + z\^ : z G (L^Z)'^}. Then the torus can be represented by the lattice cube T^v = 
{x G Z'^: \x\^ < ^{L^ — 1)} of side , equipped with the metric p oi p^. Gradient Gaussian 
fields can be easily defined as discrete gradients of Gaussian fields. However it turns out to 
be inconvenient to work directly with the space of discrete gradient fields, since the constraint 
of being curl free (in a discrete sense) leads to a complicated bookkeeping. Instead, we use 
that discrete gradient fields are in one-to-one relation to usual fields modulo a constant. To 
eliminate this constant we use the normalisation condition that the sum of the field over the 
torus vanishes. We thus denote hy X ^ the subspace 
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The forward and backward derivatives are defined as 

(y*ipYj{x) = ip''{x - Qj) - ip''{x), r = 1,. . . ,m] j = 1,. . . ,d. 

Let A: R""""^ R™^'^ be a linear map that is symmetric with respect to the standard scalar 
product (■, ■)j^mxd on r™^^'^ and positive definite, that is, there exists a constant Cq > such 
that 

{AF, F)j,,„x. > Co\\F\\l^^, for all F e R'"^'^ with \\F\\^^^, = (F, F)]/^,,. (2.4) 
The corresponding Dirichlet form defines a scalar product on ^^tv, 

{ip, := Siip, VV^(x))R™xrf, ip,ije Xn. (2.5) 

Skipping the index A^, we consider the triplet "H^ = H = of (finite-dimensional) Hilbert 
spaces obtained by equipping the space with the norms |H|_, IHI2, and \\-\\_^_, respectively. 
Here, ||-||2 denotes the ^2-norm ||v5||2 = ■, ||v^||+ = (v^,V^)+^, and ||-||_ is the dual norm 

= sup (^,</?). (2.6) 

One easily checks that ||-||_ is again induced in a unique way by a scalar product (■,■)-• The 
linear map A defines an isometry 

A-.n+^H., ^^A^ = V*{AV^). (2.7) 

Indeed, it follows from the Lax-Milgram theorem that, for each / G H-, the equation 

{ip,v)+ = {f,v) for all ven+ (2.8) 

has a unique solution (f G "H^-. Hence ^ is a bijection from to 7i . Moreover 

II^V^II- = snp{{A(p,v) : \\v\\+ < 1} = sup{(v9,t;)+ : \\v\\+ < 1} = \\(p\\ + . (2.9) 

Thus A is an isometry from 'H+ to H-. In view of the symmetry of A it follows that 

{ip,ip)- = (^"V,'^"V)+ = {A-^^,AA-^i)) = {A-^ip,ip). (2.10) 

In the more abstract construction of Brydges and Talarczyk |4] it is important that the 
operator A can be written as 

yi = s*s, (2.11) 

where S* denotes the dual of !B. This is indeed possible in our case. Since the operator A is 
symmetric and positive definite it has a positive square root A^^"^ and we can define S by 

(!B^)(x) = {A'/'V^){x). (2.12) 

This yields 

{^,i;U = {'B^,3^l;), 11^.^ = ||S^||2. (2.13) 

In the following, however, we will not use the operator S explicitly and will, instead, directly 
use that {Alp, tp) = {AVip, Vtp) and exploit that the right hand side is sufficiently local in ip 
and ip. We do, however, make crucial use of the assumption that A is positive definite in the 
proof of Lemma 13. 4[ For many of the other estimates it would be sufficient to assume that the 
operator A is positive which is implied by the weaker condition that A is positive definite on 
matrices of rank one, i.e., {A{a 6), a 6) > Co|ap|6p for all a G R"*, b G R'^. 
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Consider now the inverse Qa = -A"^ of the operator A (or the Green function) and the 
corresponding bihnear form on Xn defined by 

GAif, ^) = {Ga^P, i^) = i^)-, e Xn- (2.14) 

Given that the operator A and its inverse commutes with translations on Tat, there exists a 
unique kernel Ca such that 

{Qa^){x)= Y,CA{x-y)^{y), (2.15) 

(see Lemma 1331 below) . We write Ca G A^at, using AIat (in analogy with X n) to denote 
the space of all matrix-valued maps on T^r with zero mean. Notice that if the kernel Ca is 
constant, Ca{x) = C for any x E In, where C is a linear operator on R"*, then {Ca<p){x) = 
X^yeiAT ^a{x — y)(p{y) = C J2yeTr^ '^{y) — V ^ '^N- It is easy to see that the function 

GA,y{-) = Ca{- — y) is the unique solution GA,y G M-n of the equation 

AGA,y={6y-^,)l, (2.16) 
where 1 is the unit m x m matrix. Notice that for any a G R"^ one has: 

{AGA,y)a = {6y - ■^)a e Xn- 

We now state our main results. 

Theorem 2.1. The operator Ga'- "H- — admits a finite range decomposition, i.e., there 
exist translation invariant positive- definite operators 

eA,k: ^ {eA,k¥^){x)= J2<^AA^-yMy)^ k = i,...,N + i, (2.17) 

such that 

N+l 

Ga = ^ QA,k, (2.18) 

k=l 

and for each associated kernel Cam ^ -M.^ there exists a constant matrix Ca^u such that 

CA,k{x -y) = CA,k whenever poo{x, y) > for k = 1, . . . , N. (2.19) 

Moreover, for any multiindex a, there exist constants Ca{d) > and ri{a,d), depending only 
on the dimension d, such that 

||V"CA,fc(x)|| < C«(d)L-('^^^)('^^2+|a|)^r,(M) (2.20) 

for all X E In (ind k = 1, . . . , N + 1. Here, V" = HiLi and V° = id, and ||-|| denotes the 
operator norm. 

Note that since any function in .Ytv has mean zero the kernel CA,k = CA,k — CA,k generates the 
same operator C^^^. Thus f l2.19p indeed guarantees that C^ ^ has finite range. See Lemmas 13.61 
and 13.51 for further details. 



The operator A, its inverse C^, and the finite range decomposition itself, depend on the linear 
map A: r™^'^ — s. r'tixc! q^^, major result is that the finite range decomposition can be defined 
in such a way that the maps A i— )■ CA,k are real-analytic, as long as A is positive definite. 
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Let £sym(R"''"^) denote the space of linear maps A: pj^^^ — > -pj^xd. ^^j^g^^ qj-q symmetric with 
respect to the standard scalar product on r™^^'^ and let 

U:= {Ae rsymlR"""^) : {AF, F)^,^., > for all F G R™^^ F ^ O} (2.21) 

denote the open subset of positive definite symmetric maps. 

Theorem 2.2. Let d > 2 and let a be a multiindex. There exist constants Ca{d) and r]{a,d) 
with the following properties. For each integer N > 1, each k = 1, . . . , + 1 and each odd 
integer L > 16 there exist real-analytic maps A i— )■ CA,k from U to M-n such that the following 
three assertions hold. 

(i) // QA,k denotes the translation invariant operator on induced by CA,k then 

N+l 



Ga — ^ ^A,k- (2.22) 



k=l 



(ii) There exist constant m x m matrices CA,k such that 



CA,k{^) = CA,k z/poo(x,0) > (2.23) 



(iii) // {AqF, F)R,nxd > co||F||^„xd for all F e R"^""^ and cq > then 

j 

sup 



(V"D^C^„,fc(x)(i,...,i) <Cm(-) j!L-('^-i)('^-2+H)^.(M). (2.24) 

||A||<1 V'^O/ 

for all X E In o,nd all j > 0. Here V" = Y[i=i ' '^^ '^^^ ll^ll denote the operator 
norm of a linear mapping A: r™-^'^ — ). R"ix<i^ and the jth derivative with respect to A in 
the direction A is taken at Aq. 

The condition L > 16 can be dropped. However, in applications to renormalization group 
arguments, one needs to choose L large also for other reasons. Given 6 G (0, 1/2), the property 
(ii) can be strengthened to CA,k{x) = CA,k if Poo{x, 0) > 6L^, provided that L is large enough. 
Then the constants Ca{d) and r]{a,d) depend also on 6. We refer to Remark 13.101 for further 
details. 

The proof of the existence of the finite range decomposition, i.e., (12.180 and (12.190 of Theo- 
rem [2]T1 is given in Section |3l The remaining proof of the regularity bounds is given in Section |4] 
for a fixed A. Real- analytic dependence on A and the bounds f l2.24p are established in Section [5l 

The research of R.K. was partially supported by the grants GACR 201-09-1931, 201/12/2613, 
and MSM 0021620845. 



3. Construction of the finite range decomposition 

In this section we prove the existence part of Theorem 12.11 via an extension and adaption of 
the methods in [4] to our case. The proof of the estimates is given in Sections |4] and [H 

The existence of a finite range decomposition is contained in Proposition 13 . 81 and Proposition [3l9] 
below and their proofs are built on the following auxiliary results. 

First, for the construction of the decomposition we consider the discrete cube 

g = {i,...,/-i}'^ (3.1) 
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for some / G N, / > 3. We can identify Q with a subset of T^v once / — 1 < . Similarly, any 
shift Q + X d Ttv- For any x e Tat, consider the subspace 

-HiQ + x) = {(/? G ^: v? = in Tn\{Q + x)}. (3.2) 

We write Hj^^Q+x) and l-L^^Q+x) for the same space equipped with the scalar products (■,■)+ 
and (■,■)- 5 respectively. We denote by 11 ^ the (■, ■)+-orthogonal projection '^+{Q + 

and set Pj, = id — U^. Thus U^ip G H+iQ + x) and 

(i7,(/p,^)+ = for all ^P^n+iQ + x). (3.3) 

For any set M C Ajy, we define its closure by 

M = {x G Atv: distoo(x,M) < 1}, distoo(x, M) := min{poo(2:, y) : y G M}. (3.4) 

In particular, 

Q = {0,...,lY. (3.5) 

We also define 

Q.:={0,l,...,l-lY. (3.6) 
Lemma 3.1. For any if G we have 

(i) Ai^Pxi^) = const in Q + X, 

(ii) P^ip = ip im:iy\{Q + x), 

(iii) n^Lf = (plq+x if<^ = on {Q + x)\{Q + x). 

Remark 3.2. This shows that P^ip is essentially the A-harmonic extension in Q + x. Thus we 
would expect that P^ is (locally) smoothing and suppresses locally high frequency oscillations, 
while Tlx suppresses locally low frequencies. This will he made precise in Lemma \4-l\ below 
where we show the corresponding estimates for the averaged operators 7 = XIxst ^'^^ 
Jl = id-7. 



Proof, (i): By (13. 3p we have for all ip G 9^+(Q + x) the relation {Px'p,'ip)+ = and hence 
{A{Px(p), ip) = 0. Taking ip = 6^ — 6z for any pair of points v,z E Q + x we get A{Px(p){v) = 
A{Px(p){z). This proves (i). 

(ii) : This follows from the fact that n^ip belongs to Hj^{Q + x) and hence vanishes outside 
Q + x. 

(iii) : It suffices to consider the case x = and we write II for IIq. Let (p = iplq. Then 
(f G 1-L+{Q) and hence IICp = (p. Moreover ^ — (p vanishes in Q. Thus V(v2 — (p) vanishes 
in Q-. Hence {ip — <p,ip)+ = for all ip G 'H+{Q) since Vip is supported in Q_. Therefore 
n{ip — (fi) =0 which yields the assertion. 

□ 

Lemma 3.3. 

(i) IlxIIy = whenever {Q^ + x) fl (Q_ + y) = 0, 

(ii) IIx^ = whenever suppip fl (Q + x) = 0. 



Proof. (i): For any (p,ip & the functions il^v^ and Ilyip vanish on Tn\ {Q + x) and 

T7v\(Q+2/), respectively. Hence, V n^tp and VIIyLp vanish on Tj^\{Q_+x) and on Tjv\((5-+2/), 
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respectively. Assuming now that Q- + x and Q_ + y are disjoint and taking into account fl2.5p 
we get 

n^ny^)+ = {n,^, ny^)+ = ^ (A(vi7,^)(z), {vny^){z))^rr.., = o. (3.7) 

(ii): For ip G H+iQ + x) we have Atp = in Tat \ (Q + x). Thus for any ip G 'H+ with 
suppv? n {Q + x) = we get = {ip,Aip) = 0. In view of (13. 3p this yields Il^ip = 0. □ 

Next, consider the symmetric operator 

T = ^ E (3-8) 

on The following result is the key estimate for the finite range decomposition construction. 
Our proof is a slight modification of the argument in [1]. 

Lemma 3.4. For any ip G "H^ we have 

(i) < (n^ip, if)+< {iQ^+,AVip, v</^) , 

(ii) < {7ip,ip)+ < {(p,(p)^ and the inequalities are strict ifipj^O, 

(iii) (T<^,T(p)+<(T<^,<^)+ . 

Proof, (i): We have {nxip,ip)+ = {ip^IIx^)^ = {IIx^, n^ip)^ > 0. For the other inequality 
we use that VTI^^ is supported in Q- + x. Thus 

{n^if, ip)^ = {Avn,^, v<^) = {Avn,^, Iq^+.v^) . (3.9) 

Since A is symmetric and positive definite the expression (F, G)a '■= {AF, G) is a scalar product 
on functions Z*^ — )■ r"*^*^. Thus the Cauchy-Schwarz inequality yields 

{Avn^v, 1q_+,.V(p) < {Avn^v, vn,^y/^{AiQ_+,v^, 1q„+,v<^)^/2 

= {n,^,n,ip)'f{iQ^+,AVip,Vipy/\ (3.10) 

Together with (13.91) this yields the assertion since {nx(p,(p)+ = {Ilxip, nx(p)+. 

(ii) : Since J^xeT^ ^Q-+^(y) ~ ^'^ all ?/ G Tat the inequalities follow by summing (i) over 
X G Tat. If (T(y9, ip)j^ = then {Ux^p, (p)+ = for all x G Tat and thus n^^p = and Px^p = 
Lemma [371] implies that there exist constants such that {Aip){y) = for ally E Q + x. Since 
/ > 3 the cubes Q + x and Q + {x + ei) overlap and this yields Cx = Cx+ei ior all i = 1, . . . ,d. 
Thus Cx is independent of x. Since Aip E this implies c = 0. Hence A(p = and therefore 
ip = 0. 

Now suppose that {7ip,ip)+ = {(p,(p)^. This implies that for all x G Tat we have {nx<p,<p)+ = 
{^Q-+xA'Vip,'Vip) . We claim that the last identity implies that Vip{x) = 0. Indeed, if 
1q_+xVv7 = we are done. Otherwise the identity can only hold if the inequality in (I3.10p is 
an identity. In particular we must have Vil^v? = ^l-g.+x^f and A = 1. Now il^v^ vanishes 
outside Q + X and in particular at the points x and x + e^. Thus VUx^p^x) = and hence 
Vip{x) = 0. It follows that ip is constant on Tat and hence ip = since ip has mean zero. 

(iii) : It follows from (ii) that {(p,ip)^ := (7(p,ip)^ defines a scalar product on Thus the 
Cauchy Schwarz inequality and (ii) yield 

(Ty;,^)+ < (T^,(^)f (T7/;,^)f < (Ty.,y.)f (7/;,^)f . (3.11) 
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Taking ip = 7(f we obtain the desired estimate. □ 
Consider the operator 7' : — ?■ 'H_ dual with respect to T and defined by 

(tV,^/') = (<^,TV;), ip en-,i/j en+. (3.12) 

Notice that 

7' = A7A-\ (TV,^)- = (<^,TV)-, and (TV,<^)- = (Tyi~V,'^~V)+- (3-13) 
Indeed, for any ip G 'H+, we have 

{7'Ap,,ij) = {Ap,,7ij) = (<^,TV;)+ = {7^,^P)+ = {A7^,^P), (3.14) 
and this yields the first identity in fl3.13p . Now 

(TV,V^)_ = {A-^A7A^^^,tlj) = (Tyi-V,^^~V) = (T^"V,'^"V)+- (3.15) 

Since the last expression is symmetric in if and ip we get the second identity in (13.131) and 
taking ip = ip we obtain the third identity. Similarly, we have U'^ = ATI^A'^ for the dual of 
n^. Notice that 

n'^ip = whenever suppc/? H {Q + x) = 0. (3.16) 

Indeed, considering any test function ip E Xjy, we have {n'^ip,tp) = {ip,n.j.ip) = 0. We also 
consider the operator 

Ji:=id-7 and its dual 3^' = id - T'. (3.17) 
It follows from Lemma [3.4( ii) and (I3.13P that 

(TV,</^)->0, (3lV,V')->0, (a'V,TV)_ < (^V,^)- for all G \ {0}. (3.18) 

We next discuss the locality properties of translation invariant bilinear forms, operators and 
the corresponding kernels. These properties would be obvious if we consider bilinear forms on 
Vat since then we can use the Dirac masses 6x as test functions. Dirac masses, however, do 
not belong to Xn and hence we need to use test functions with broader support which makes 
the conclusion of Lemma 13.61 below nontrivial. We begin by recalling the relation between 
bilinear forms, operators and kernels. The translation operator Ta,a G Tat, is defined by 
Ta^piu) = (p{y — a), so that TaSy = Sa+y. Recall that AIat denotes the space of all matrix-valued, 
periodic maps with zero mean. 

Lemma 3.5. Let B be a translation invariant bilinear form on Xj^, i.e., 

Bijaip^Taip) = B{ip,ip) for all G AT AT, for all a G T^r. (3.19) 

Then the following assertions hold. 

(i) There exists a unique linear operator S: X n X n such that 

(!B(^, i)) = B{p, ij) for all ip,tp e X^. (3.20) 

Moreover S is translation invariant, i.e., "Bra = Ta'B. 

(ii) There exists a unique matrix-valued kernel B G M-n such that 

{'Bif){x) = ^B{x-y)(f{y) for all x G T^r, for all (feX^. (3.21) 
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Moreover for i5: Tjv — )■ r™-^"* we have 

= ^ - y)<^{y) for all x G for all <^ G Xn- (3.22) 

if and only if 

B-B = C (3.23) 

with a constant m x m matrix C . 

(iii) If B' E X N denotes the kernel of the dual operator 'B' then 

B\z) = B{~z). (3.24) 

(iv) //Si and S2 are translation invariant operators on X and 'B^ = then S3 is 
translation invariant and the corresponding kernels Bi G Ad^ji = 1,2,3, are related by 
discrete convolution, i.e., 

Bs{x) = {Bi * B2){x) := ^ Si(x - y)B2{y) for all x G T^v. (3.25) 

y6Tiv 

Proof. We include the elementary proof for the convenience of the reader. 

(i) : Existence and uniqueness of S follows from the Riesz representation theorem. To prove 
translation invariance let D := Sr^ — TqS. We have 

CDif, Ta^) = {BTaif, Ta^) - (r,S(/^, Ta^) . (3.26) 

Since is an isometry with respect to the scalar product (■, ■) we get 

CDip, Ta^) = B{Ta^, Ta^j) - {B^, ^j) = B{TaV, T,^) - B{^, i,) = 0. (3.27) 

This holds for all ip^ip G X n. Hence D = 0. 

(ii) : To show the existence of B note that (5o — ^ro)*^ ^ X ^ with any a G R'" and define 

S(x)a:=S[a((5o-;^)](x). (3.28) 

Using that rj,(5o) = 5^,, we get 

B{x - y)a = S[a(5o - ■j)ni)]{x - y) = ry(S[a((5o - ;^)])(x) = 

= Sr,[a(5o - i^)](a;) = S[a(5, - x^)](x). (3.29) 
Observing, further, that for any ip E X ^ we have 

we get 

(a5^)(x) = nv{y){5y - jh)]{x) =Y.B{x~ y)^{y). (3.31) 



This shows the existence of B. Suppose now that Bipi^x) = XlyeTjv ~ y)^iy) aA\ ip E X 



N 



and set C{z) = B{z) — B{z). The choice (p = {6^ — So)a with an arbitrary a G R™" yields 

C{x - ei) = C{x) for all x G Tat, for all z = 1, . . . , d. (3.32) 

Thus B — B = C with a constant m x m matrix C. If in addition B G AIat then this implies 
B = B. Conversely if i3 = B+C then B and B generate the same operator since X^xgtat fi^) — 
for p G ATjv. 
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fiii): We have 



= ^ (S(-[x-2/])<^(2/),^(x))r™. (3.33) 

Hence = Xlyei '^("[^"^/Dv'l?/) the uniqueness result in (ii) imphes that B\z) = 

B{-z). 

(iv): One easily verifies that Si *i32 G M.n and that {'Bi'B>2^){x) = J^yeT^ * B2{x — y)ip{y). 
Thus the assertion follows from the uniqueness result in (ii). □ 

For two sets Mi, M2 C Tat we define 

distoo(Mi,M2) := mm{p^{x,y): X E M^y e M2}. (3.34) 

Lemma 3.6. Let B be a translation invariant bilinear form on X ^ and let 25 and B G A^n be 
the associated operator and the associated kernel, respectively. Let n be an integer and suppose 
that > 2n + 3. Then the following three statements are equivalent. 

(i) B{ip,ip) = whenever distoo{suppip,suppip) > n. 

(ii) There exists an m x m matrix C such that B{z) = C whenever poo{z, 0) > n. 

(iii) supp!B(y9 C supp(/9 + {—n, . . . , n}'^ for all ip E X^- 

Proof. The implication (ii) =^ (iii) is easy. Set B{z) = B{z) — C. Then B{z) = if 
Poo(^) > n with Poo{z) = Poo{z,0) and by Lemma [375l (ii) we have 

(a3^)(x)= Y.B{x-yMy). (3.35) 

J/GTjv 

If X ^ suppv9 + {— n, . . . , n}'^ then either y ^ suppip or y E suppy^ and Poo{x — y,0) > n. In 
either case "Bipi^x) = 0. 

The imphcation (iii) =^ (i) is also easy. Suppose that distoo(suppv9, suppV') > n. Then (iii) 
implies that distoo(suppS(y9, suppz/^) > 0, i.e, "Bip and ip have disjoint support. Thus B{ip,ip) = 
(!B<^,^) = 0. 

To prove the implication (i) =^ (ii), consider the torus Tat with the fundamental domain 
Ljv = {-^^, • • • , and set 

M := {-n,...,ny, (3.36) 
M_ := {— n, . . . , n + l}'^ , and the closure of M, (3.37) 
M = {-(n + l),-n,...,n + l}"*. (3.38) 

Note that by the assumption > 2ri + 3 the set T^r \ M is nonempty. 
We first show that for alH, j G {1, . . . , d} we have 

ViV*B = in \ M. (3.39) 
To see this let ^ G Tat \ M and consider 

= (<5«+e, - <5s)a, ^ = (5e, - So)b with a, b E R™. (3.40) 
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Since \ei — ej\oo < 1 we have 

distoo(suppv9,suppz^) > PooiO,^) -1 >n + l. (3.41) 

Hence 

= Bicp, ^) = ( + e, - e,) - + e,) - - e,) - b) = (V,,V*S(Oa, b) (3.42) 
for every a,b E R™ and thus 

V,V*S(0 = 0. (3.43) 
Next we show that for j G {1, . . . ,d} there exist a matrix Cj such that 

V*B = Cj in \ M. (3.44) 

Fix J and set / = V*B. Using the shorthand I = {—n — 1, . . . , n + 1}, let 

XI G {-^, . . . , ^} \ /, x' := (X2, . . . , X,) G {-i^, . . . , ^j-^-^ (3.45) 

Then x = {xi, x') G Tat \ M and hence for z 7^ 1 we have (Vi/)(xi, x') = 0. Thus there exists a 
matrix- valued function gi on { — . . . , ^ ^"-^ j \ J such that 

/(x) = (7i(xi) ifxi^J. (3.46) 

In the same manner, there exists a function g2 such that f{x) = (72 (X2) if X2 ^ /. This eventually 
implies (K^ . 

Further, 

V*fi(-n - l,x') - V*i3(-n - 2,x') = {ViV*B){-n - 2,x') = 0. (3.47) 

Hence we also have V*B{x) = Cj if xi = —n — 1. Arguing similarly for the other components 
of X we get 

W*B = Cj inTjv\M_. (3.48) 

We now show that Cj = 0. Assume without loss of generality j ^ 1. For x\ ^ {l \ {— — 1}) 
we have 

V*i5(xi, x') = Cj for all x' G {-^, . . . , ^}'^-^ (3.49) 
On the other hand X]^,^^ ^^-1 V*i3(xi,x') = since B is periodic and j 7^ 1. Thus 

C,=0. 

Arguing as in the derivation of (13.441) we conclude from (13.481) and the fact that Cj = that 
there exists a matrix C such that 

B = C in Ttv \ M_. (3.50) 

In addition we have 

B{n + 1, x') - B{n + 2, x) = {VlB){n + 2, x) = 0. (3.51) 

Hence B{x) = C for xi = n + 1. Arguing similarly for the other components of x we get B = C 
in Tjsf \ M. This finishes the proof of Lemma [3.61 (ii). □ 

Lemma 3.7. Suppose that distoo(suppv9, supp?/;) > / — 1. Then 

{Tip, tp) = 0, (TV, V^) = 0, {Jiifi, tp) = 0, (3^V, ^) = 0. (3.52) 
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Proof. It suffices to prove the first identity. The second follows by exchanging ip and if) and 
the third and fourth follow since = id — T and 3^' = id — T'. By Lemma [3.31 we have 

n^ip = if supp(y9 n (Q + a;) = (3.53) 

and it follows from the definition of 77^ that suppil^v^ <Z Q + x. Assume (Ty?, ip) ^ 0. Then 
there exist x G T^r such that {Uxif, ip) ^ 0. Thus supp'?/'n((5+x) 7^ and suppv9n((5 + a;) 7^ 0- 
Therefore there exist ^ E Q and ( E Q such that x + G supp(/?, x + C G supp?/^. Thus 

x + ^-{x + C) = ^-(e{-{l-l),...J-lY. (3.54) 

Hence distoo(suppv9, supp-?/') < / — 1. □ 

Consider now the inverse C = A''^. (For time being, we omit the reference to the matrix A in 
the notation for C. We will reinstate it in Section [5|, where we explicitly discuss the smoothness 
with respect to A.) The main step toward the decomposition, is to subtract a positive definite 
operator from C in such a way that the remnant is positive definite and of finite range. We 
define 

Gi := e - JieJl' , which yields Ci = C - 7^ * C * 7^'. (3.55) 
Proposition 3.8. Both Ci and JlGJl' are positive definite and Ci has finite range, i.e., 

{Qiip,ip) = z/ distoo(suppv9, supp?/;) > 2/ — 3. (3.56) 
In particular, there exists an m x m matrix C such that 

Ci{z) =C if p^{z, 0) > 2/ - 3. (3.57) 



Proof. For any (p,ip G Xn, we use f l2.14p to get 

{JlQJi'ip, ip) = {Ji'ifi, > 0. (3.58) 

If Jl'ip = then fl3.18p implies that = 0. Thus 'RQ'Jl' is positive definite. Furthermore, 
{Grip, ij) = {Qif, ij) - {GJl'ip, Jl'ij) 

= (<^, ^A)_ - (Ji'ip, 3^V)_ = (TV, + (<^, TV)_ - (TV, T»_. 
Thus (13.181) implies that Ci is positive definite. 

To evaluate the range of the quadratic form (61^9, ip), we inspect the terms on the right hand 
side of f l3.59p . For the ffist (and similarly the second) term, we have 

(TV, = ^ E (^^^' = ^ E (^^^' ^^^)- (3.60) 

In view of fl3.16p . a term in the sum vanishes at x except when the supports of (p and ip both 
intersect Q + x. Therefore, the scalar product is zero whenever the distance of the supports is 
strictly greater than / — 1. The second term of the bilinear form Gi{(f,ip) := {Ci(p,ip) is the 
double sum 

(^V,T»_ = ^ E ^ E (3-61) 

By Lemma 13731 we have U'^Uy = ATIxHyA'^ = whenever (Q- + x) fl {Q- + y) = 0, i.e., if 
Poo{x,y) > / — 1. Hence the double sum only contains a non-zero contribution if there exist x 
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and y such that poo[x,y) < / — 1, suppy? fl Q + x 7^ 0, and supp?/' HQ + y ^ 0. Hence there 
must exist ^ Q such that x + E suppip and y + ( E suppip. Hence 

distoo(suppv9,suppz/;) < Poo{x + ^-{y + C),0) < poo{x-y,0)+poo{^-C,0) < Z-1 + /-2 < 2/-3. 

(3.62) 

This proves fl3.56p . and (13.571) follows from Lemma [3.61 □ 

We construct a finite range decomposition by an iterated application of Proposition 13. 8[ Let 
L > 16 and consider 

Qj = {l,...Jj- lY with Ij = [|L^J +1 for J = 1, . . . , AT. (3.63) 

Here [aj denotes the integer part of a, the largest integer not greater than a. In particular we 
have 

IU<1,<IU + 1. (3.64) 
We define Tj, 7'^ and Jl'j as before with Q replaced by Qj and set 

Cfc := (3^1 . . . Jik~i)G{'R'i^^i ■ ■ ■ 3^1) — (3^1 . . . ^k-i^k)G{'Ji'f^Ji'i^_i . . . ^R'x), k = 1, . . . , N, (3.65) 
and 

Ca^+I '■= (^1 • • • ^N~l ■ ■ ■ 3iAr)C(3^7V"'^Ar-l • • • ^l)- (3.66) 

With these definitions, we show that the sequence {Gk}k=i,...,N+i yields a finite range decom- 
position. 

Proposition 3.9. Suppose that L > 16. Then the operators Qk satisfy 

(i) e = Ef=Y e.. 

(ii) Cfc is positive definite for k = 1, . . . , N + 1. 

(iii) For k = 1, . . . , N the range of is bounded by , i.e., 

{Gk^,ip) = z/ distoo(supp(/9, supp?/;) > ^L'' (3.67) 
and there exist m x m matrices Ck such that 

Ck{z)=Ck %f Poo(^,0) > \LK (3.68) 

Remark 3.10. 

(i) Let 5 G (0, 1/2). Then we can obtain the sharper conclusion Ck{z) = Ck if p{z,0) > 5L^ , 
provided that L is large enough and we choose the integers Ij sufficiently small, e.g. we may 
take Ij = [SL^/4:\ +1 if L >8/6. Of course the regularity estimates (12.201) then depend on 6 
and degenerate for 5—7-0. 

(ii) The restriction L > 16 can be removed. If 6 < L < 15 we can take li = 3 and for j > 2 
define Ij as before. One can easily check that in this case we still have —1+2 X]j=i(^i^l) — L''/2 
and this yields the desired assertion (see the proof below). If 3 < L < 5 one can skip the first 
few renormalization steps. Formally one can take li = I2 = 2 and define Ij as before for j > 3. 
Then 7^=72 = 0, O^i = = id, Ci = 62 = and -1 + 2 Ei=3(^j - 1) < ^72- 

Proof. Assertion (i) follows directly from the definition. To prove (ii), set 

■■= K-i ■ ■ ■ 'ipk ■■= K-i ■ ■ ■ ^1^, k=l,...,N + 1. (3.69) 
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Inductive application of f l3.18p shows that (fk = imphes (f = 0. Now, directly from definitions, 
'■p) = (v^7v+i, V^AT+i)-- Thus Qn+1 is positive definite. For /c = 1, . . . , iV we have 

(efc<^, ^) = ((e - 'JikQK)vk, ^k). (3.70) 

Hence by Proposition 13.81 we get (Cfc(y9, if) > Q with equality only holding if tpk = 0, which 
implies (p = 0. Thus is positive definite. 

(iii): In view of the equation {Qk'P, ^) = ((C — ^kG'Jl'i^)ipk, ipk), Proposition 13.81 implies that 

(CfcV?, i>) = if distoo(suppv2fc, supp^^fc) > 24 - 3. (3.71) 
Iterative application of Lemma 13.71 and Lemma 13.61 yields 

fc-i 

suppypfc C suppv? + {-Uk, . . . , UkY, snppipk C suppV^ + {-Uk, UkY, Uk = ^{Ij - !)• 

i=i 

(3.72) 

Thus 

k 

{Gk^p, ijj) = if distoo(suppv9fc, supp^fc) > -1 + 2 ^(/j - 1). (3.73) 

i=i 

Now since Ij - I < and En=o^~" < 2 we get ^Yj^^Sj "1) < ■ This finishes the 
proof. □ 

4. Estimates for fixed A 

To prove the regularity bounds of Theorem 12.11 and Theorem 12.21 we derive estimates for the 
Fourier multipliers of the relevant operators. To this end, we first extend the space Vat to the set 

V N = (C™j of complex- valued vectors with the subspace defined, again, as the subset 
of functions if G Vn with vanishing sum, Y1ix(^tm V^^) ~ ^- Various scalar products and norms 
are extended to complex-valued functions in the usual way, {(fiyip) := X]x-eT]v('/'(^)' '^*(?/))c'"5 
(p,il^ & Vtv, with ip*{x) denoting the complex conjugate of ip{x). 

To introduce the discrete Fourier transform, consider the set of (scalar) functions fp{x) = 
Q^{p,^)^p G Tat , labelled by the dual torus 

?^ = {p=(pi,...,p.):p,G{i^^,i^^,...,0,...,(^^ (4.1) 

The family of functions IL"'^/^}^^- forms an orthonormal basis of C^'^'' . For any ip G Vn, 
we can define the Fourier transform component-wise, 

Hp) ■= Yl P ^ (4.2) 

For ip e AT TV, we get ip{0) = J2xe7„ "ipix) = with the inverse 

^(a:) = L-^^ J2 fpi^Mp) (4.3) 
peTiv\{o} 

and 

(V9,^)=L-^'^ J2 {m,kp))c- (4.4) 

p6Tjv\{0} 
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In the same way (component-wise) we get the Fourier transform for a matrix-valued kernel 
/C e Mn, 

JC{p) := Yl fpi-^M^)- (4-5) 

xGTjv 

For a translation invariant operator %: X ^ X ]^ with the kernel /C G A^at, we get 

XV^(p) = IC{p)i^{p). (4.6) 

Indeed, 

^(p) = E E - y)mfp{-^) = E E - y)^iy)fpi-i^ - y))U-y) = 
= E E ^iy)fpi-y) = (4.7) 

Henceforth, we call /C the Fourier multiplier of X. Applying the equality fl4.6p with = afp, 
p 7^ 0, we get 

Xafp = jC{p)afp. (4.8) 

Indeed, taking into account that fp{p') = L^'^6p^p>, we have %afp{p') = KL{p')afp{p') = 
L^'^5p^p>Z{p')a = L^Hp^p>Z{p)a and thus 

%afp{x) = L-^'' Yl fp'i^)^pip') = E fp'i^)^,p'ma = Hp)afp. (4.9) 

p'6Tjv\{0} p'eTjv\{0} 



Notice also that, by Lemma [3.51 the kernel of a product of two translation invariant operators 
is given by the discrete convolution of the kernels and thus 

K:^2{p) = K.M^2{p)- (4.10) 

Now, we will study the Fourier multipliers of the operators A = V*A'V as well as the operators 
T and H introduced in the previous section. Given that Vj/p = qj{p)fp with qj{p) = e*^-'' — 1 
and V*/p = q*ip)fp, j = 1, . . . , we have 

(/la/p)r = E^i(^')^^J>.'="^^fc(^')/p for any a e R"", (4.11) 

where Arj-s,k are matrix elements of A, (A{a'^q))rj = J2T=i Ylk=i ^T,j;s,k0.sqk- In view of (14. 8p . 
we get, for the Fourier multiplier A{p) of the operator A, the expession 

{A{p)a)r = Y^W^r,j,s,kasqk{p) for any a G R™. (4.12) 

Alternatively, we can express this in terms of corresponding quadratic form, 

{A{p)a,h)c^.^ = {A{a® q{p)),h® q{p))cmxd for any a, 6 G R"". (4.13) 

It follows that the multiplier A{p) is Hermitian and positive definite. More precisely, for any 
p G Tjv \ {0}, we have 

||^(p)|| < \\A\\ \p\' and WApY'W < (4.14) 
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Indeed, the first bound follows directly from f l4.13p and the definition of the operator norm \\A\\ 
of the linear map A, 

\\A\\ := max{\AF\ : F G 0"*^^ |F| < 1}. (4.15) 

We also took into account that < |pp as follows from the upper bound in the estimate 

< |e** — Ip < valid for all t G [— 7r,7r]. To get the second inequality, we use the lower 
bound (12.41) as well as the lower bound above, to get the estimate {A{p)a,a)c"^ > ^Co|pp|ap 
for any p E In- 

For the Fourier multiplier of T defined in (13.81) . we first recall that the translation operator 
Tx is defined by {Txip){y) = ip{y — x). Hence we have 11 ^ip = TxnQ{T_xip). Note also that 
/p = e'^^'P^fp. Writing 77 shorthand for iJo, we get 

f{p)a = l-'n{^p){p) (4.16) 
for any a G C"*. Indeed, applying 7 to afp, we get 

I'n^^UM = Yl n^afpM = Yl n{T..,afp){y-x) 

= Y n{e'^^^^^afp){y - x) = ^ n{afp){y - x) e'^'^-y^ e'^^'y) 

= Yl n{afp){z)e-^^P^^^ S^^y^ = ff{^p){p)fp{y). (4.17) 

Thus, 7{afp) = l~'^n{afp){p)fp implying the claim by comparing with (14. 8p . 

We now use the symmetry and boundedness of T, with respect to the scalar product (■, ■)+' 
to deduce the corresponding properties for T(p). According to (14.81) . we have 

(7(a/p), bfp)+ = {A7{afp), bfp) = L'"'{A{p)f{p)a, b)c^ (4.18) 

for any a,b E C". Combining this with the fact that T is Hermitian with respect to (■, ■)+ and 
with Lemma [3.4( ii). we infer that 

{A{p)f{p)a,b)c-^ = {a,A{p)f{p)b)cm and < {A{p)f{p)a,a)cn^ < {A{p)a,a)cm (4.19) 

for alia, 6 G C". Since A{p) is Hermitian and positive definite, it has a unique Hermitian 
positive definite square root Ai^pY^"^ with inverse A{p)-^^^. Applying (Km . we get 

\\A{pY/'\\ < \\A\\'/'\p\ and \\Aip)-'/'\\ < (4.20) 

Setting, finally 

f{p) := A{pY^'f{p)A{pr'^' = A{pr'/'{A{p)f{p))Aip)-'/' (4.21) 

and 

jZ(p) := i:(p)V2^(p)l(p)-V2 _ J _ f (p)^ (4 22) 

we get the following bounds. 

Lemma 4.1. The operators T{p) andTZ{p) are Hermitian (with respect to the standard scalar 
product on C^) and satisfy Jor any p E Tn\ {0}, the following bounds: 
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(i) there is a constant c < oo (which depends only on \\A\\ and cq, and the dimension d) 
such that 

\\l-np)\\ = ||f(p)|| <min(l,c(b|0'), (4.23) 

(ii) there is a constant c < oo (which depends only on \\A\\, cq, and the dimension d) such 
that 

||^(p)||<mm(l,f(^ + l)), (4.24) 

These estimates show that 7 suppresses low frequencies, while 51 suppresses high frequencies, 
reflecting the idea that P^, = id — II^ is a (locally) smoothing operator (cf. Remark | 



Proof. From fl4.19p . we get 

{T{p)a,b)cm = {a,T{p)b)cm and < {T{j>)a,a)cm < {a,a)cm. (4.25) 

for any a,b G C". The operators T{p) and TZ{p) are thus Hermitian and < T{p) < 1 and, 
equivalently, < TZ{p) < 1. This implies that 

||f(p)||<l, ||^(p)||<l. (4.26) 

(i): In view of fl4.26p we may assume that \p\l < 1. We first estimate the norm of the Hermitian 
matrix A{p)T{p). First, we show that 

l''{Aip)fip)a,a)c^ = \\n{aU)\\l (4-27) 

To see this, we start from the right hand side, 

||i7(a^)||^ = {n{af,),af,)+ = {n{af,),A{af,)) = L"^'^ (77(^5), ^(5^)) = 
= L-^-^ Yl {n{^p){p'),Ap)al{p'))cr. = {n{^,){p),A{p)a)c^ = l\A{p)np)a,a),^. 

(4.28) 

Here, we first used the fact that 77 is an orthogonal projection with respect to (■, ■)+, passing to 
the second hne we used the equation A{afp){p') = A{p)afp{p') obtained as the Fourier transform 
of f l4.8p for A, then the fact that fp{p') = L^'^Sp^p' and, finally, we applied the equation fl4.16p . 

Applying ( 14. 8 P and using that afp E Xn (for p 0) and thus also n{afp) E X^, we get 
\\n{afp)\\l = {n{afp),A{afp)) = {n{afp),A{p)afp) = {n{afp), A{p)a{fp - 1)). (4.29) 

Further, given that n{afp) is supported in Q and \ fp{z) — 1| < \/d\p\l for z G Q,we have 
\\n{afp)\\, = ||i7(a/p)||,^(Q) and ||(/p - l)!^!!^ = \\fp - < Vd\p\l'/'+\ With the help of 

Schwarz inequality and this observation, we get 

\\niafp)\\l = {n{afp),A{p)a{fp-l)lQ) < ||i7(a/,)||,^(Q)M(p)|||a|v^|p|/'^/2+\ (4.30) 
The Poincare inequality ( [HI HO] ) implies that 

< c/||Vi7(a/p)||,^(Q) < -^/||77(a/,)||^ (4.31) 



with a suitable constant c. Combining this with (I4.30p and (14.140 . we get 

A 



n{ajp)\\^ < Vdc^^\p\V/'+'\a\. (4.32) 



Cn 
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Given that A{p)T{p) is Hermitian, we get 

\\A{p)fip)\\ = max{A{p)f{p)a,a)c^ < dc\\\A\\yco)\p\n\ (4.33) 

|a|<l 

The assertion follows using the second inequality in (14 .Mp . 

(ii): In view of (14.261) . it suffices to consider the case \p\l > 1. Again, we first estimate 
\\A{p)n{p)\\. Since n{p) = 11 - r(p) we get from KTT^ 

l\A{p)n{p)a, a)c.n = l\A{p)a, a) - {n{afp), afp) + . (4.34) 

Let a; be a cut-off function such that 

u{z) = 1 ii z eQ\Q, w = if dist(2,g\g) > 1 + j^, < w < 1 and |Vu;| < c|p| 

(4.35) 

with a suitable constant c. By Lemma IXT] (ii) we have 77(1 — uj){afp) = (1 — uj)lQafp. Hence 
in{afp), afp)+ = {n{aufp), afp)+ + (a(l - uj)lQfp, afp)+ (4.36) 

and 

(a(l - u)lQfp,afp)+ = (a(l - u)lQfp,A{afp)) = 

= {a{l-uj)lQfp,A{p)afp)= J](l-a;) (l(p)a,a)c".. (4.37) 

zeQ- 

Here, in the last step, we used that J^zeqi^ — u) = J2zeQ-(^ ~ ^) since w = 1 on Q_ \ Q. 
Using that |Q_| = Z"', the equations (Oil) and (06]) with KST^ yield 

l'^{Aip)nip)a, a)cm = -in{aufp), afp)+ + ^ u;(l(p)a, a)c-. (4.38) 

26Q- 

Given that u is supported in a neighbourhood of order 1 + around the boundary Q\Q oi 
Q, the last term is easily estimated 

I V uj{A{p)a,a)c^\ < Adl'^^^l + l-)\\A{p)\\\a\' < ^dl^'^l + ^)\\A\\\p\'\a\'. (4.39) 
^ \P\ \P\ 



z&Q- 



To bound the remaining term we introduce another cut-off function u that satisfies the following 
conditions, 

2 3 

uiz) = 1 if dist(z, Q\Q) < 2+ — , = if dist(z, Q\Q) > 3+—, < cj < 1 and |Vw| < c\p\. 

\P\ \P\ 

(4.40) 

Then 

{n{aujfp), afp)+ = {aufp, n{afp))+ = {aufp, n{aufp)+ + {aufp, (1 - uj)lQafp) + 
= {aufp, n{aufp)+ + {A{aufp), (1 - cI;)lQa/p) 

= {aufp,n{aufp))+ (4.41) 
since 1 — a; and A{aufp) have disjoint support. Thus 

3 

\{n{aufp),afp)+\ < ||aa;/p||+||i7(a(:;/p)||+ < ||aa;/p||+ ||a(:;/p||+ < f\\A\\\p\H'^^\3 + —)\a\^, 

(4.42) 
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where we used that u and a) are supported in a strip of size 3 + 3/|p| around Q \ Q, that 
< I, and that the gradients of u, u and fp are bounded by c\p\ and the constant c in 
( I4.42P is suitably chosen in dependence on c and d. The combination of fl4.38p . fl4.39p and f l4.42p 
now yields the estimate 

\\Ap)np)\\ < (!c + 4d)||A||bpl(l + p). (4.43) 

In view of f l4.14p this finishes the proof of (ii). □ 
As in the previous section assume that L > 16 and consider 

Qj = with Ij = [|LJJ + 1 for J = 1, . . . , iV. (4.44) 

Also operators Tj, 7'^ and Jl'j, as well as C^, k = 1, . . . , N + 1, are defined as before (cf. (13.651) 
and dmi ). 

We define A^^'^ via the action of the corresponding Fourier symbol, A^/'^(p{p) = A^^'^{p)ip{p). 
Similarly we define A^^^^. Then the operators ■= A^I^'R^A"^!'^ and Tfc := yi^/2j^y^-i/2 
Hermitian. For = 1, . . . , TV define 

Mfc := ^Ri . . . ^Rfc_i^fc and Mq := id. (4.45) 

Since = id - and ^ = A'^ = A'^I^A'^I'^ we have for A; = 1, . . . iV 

Cfc = A~^l''\M.k^^jk^_^-^^M!^\A~^l'' 

= yi-i/2[(TfcM,_i)Ml_i + Mfe_i(T,Mfc_i)' - TfeM,_i(?fcMfc_i)']yL-^/2 ^4 45) 

and 

e^+i = A-^'^MnM'nA-^'^ . (4.47) 

To estimate the corresponding kernels Ck and their derivatives we use formula (14. 3 p for the 
inverse Fourier transform. This yields 

sup ||C,(x)|| < -4 E 11^^(^)11 (4.48) 

p6Tjv\{0} 

and for any multiindex a = (ai, . . . , a^), 

^sup ||V"C,(x)||<-^ J2 ||C,(p)g"|| < (4-49) 

P6Tjv\{0} p6Tjv\{0} 

where qj = e*^^ — 1 and q"" = Y['j=i I'j^ ■ 

Finally, a bound on ||Cfc(j9)|| will be based on (I4.46P (respectively, (I4.47P ) combined with 
bounds on ||A^fc(^»)|| and \\fk+iip)Mkip)\\. 

The estimate of the latter will depend on Namely, slicing the dual torus into the annuli 
A, := {pETnX {0} : nL^^ < \p\ < 7rL-^+'}, j = 1, N, (4.50) 
with the complement 

Ao:={peTN\{0}:\p\>7r}, (4.51) 
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and defining, for any c > 1 and j, k = 0,1, . . . , N , the step functions 
and 

[_ L(fe-j)(fc-J + l)/2) II t /Ij, J ^ ft, 

we liave the following estimates. 

Lemma 4.2. There exists a constant c (depending only on cq, \\A\\, and d) such that for any 
odd L > 16 and any N en,N >1, 

\\Mk{p)\\ < M,^,,l{p) for k = 0,...,N, (4.54) 

and 

\\fk+i{p)Mk{p)\\ < Mk,cAp) fork = 0,...,N-l. (4.55) 

Proof. Let p G Aj. For k < j both bounds follow from the bounds ||7^„(p)|| < 1 for 
n = 0,1,..., k and ||7fc+i(p)|| < c(L^"'"^|p|)^. Now, assume that k > j and recall that (after 
increasing the constant c from f l4.24p by a constant factor) 

ll^n(p)ll<:^(l + p)<X^' ioTn = j + l,...N. (4.56) 

Thus 

n iiK..(P)ii < n ii«..(p)ii £ n = L,.-mLm - f^-") 

n=l n=j+l n=j+l 

The first estimate for k > j follows, with the second estimate implied since HTfc+id?)!! < 1. □ 

To combine these bounds for an estimate on the Fourier multpliers Ck (p) , we use the following 
Lemma. 

Lemma 4.3. Letn be a nonegative integer andc > 1. Then there exists a constant d (depending 
on parameters c, n, and the dimension d ) such that with 

rj = rjin, d) = max(i(rf + n- if, d + n + 6) + 2 (4.58) 

and for all integers L > 3, N > 1, and all k = 1, . . . , N + 1, we have 

^ J2 M,_i,e,L(p)M,-i,c,L(p)|pr-'<c'L''L-('="^)('^+'^-2), (4.59) 
peTjv\{0} 

^ A4-i,c,L(p)Mfc-i,c,L(p)br-'<cc'L''+«L-('=-i)('^+"-2), (4.60) 

P6Tjv\{0} 

^ J2 MN,,,L{p)M^^,,Lip)\pr'<c'L^L~''id + n-2). (4.61) 

peTjv\{0} 
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Proof. It suffices to prove the ffist bound for k = 1, . . . N + 1. The second and third bounds 
follow employing the inequalities Mk-i,c,L < cL^ Mk-i,c,L and Mn^c,l < Mn^c,l, respectively. 

To prove the ffist estimate we split the sum into the sum of contributions over the annuli Aj. 
For p & Aj, we have 

< 7r"-2rf"/2L2L(-^'+i)("-2). (4.62) 

Indeed, for j 7^ 0, we get 

|p|n-2 ^ 2,"^ax((2-n),0)_^(-i+l)(n-2)^n-2 ^ ^„_2_^2_^(-j+l)(n-2) _ {4.Q3) 

The expression {2 — n) in the term max((2 — n), 0) stems from the fact that for n = and n = 1, 
we actually employ the lower bound on \p\ from (14.501) . For j = 0, we have tt < \p\ < \/d\p\oo < 
Vdir and thus Ip]""^ < Tr'^-^d"/^ < 7r"-2d"/2L". The size of the annuli can be bounded as 



< 7r'^L(-J+^)'^. (4.64) 



IdN 

As a result, for j > k — 1, 

^ 5^ M,„i,,,,(p)M,_i,e,,(p)br-^ < c7r"+'^-2d"/^L2L(-^-+i)('^+"-2)L«L^('=-i-^-) < 

^ ^^n+a!-2^n/22^8+a!+n2^-(fc-l)((i!+n-2)2^-(j-(fc-l))(d+n+2)_ (4_g5) 

and 

1 ^ ~ 



with c = 2c7r""'''^"2(i"/2 since 



^ i-0-,.-i„,.+..«, _ ^ i-«.+„«, ^ ^_J^. (4.67) 
Now consider j < k — 1. We get 

^ E M,_i,,,(p)M,_l,,,(p)br-^ < ^n+.-2^n/2^2^(-,>l)(.+.-2) ^g^-^^-^^^^ < 

< ^n+d-2^nl2-^2j^{-k+l){d+n-2)^ik-m+n-2) f ^^"'"'^ (4.68) 



Setting j' = A; — 1 — j we get 

k~2 k~\ 



(4.69) 
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Consider the integer £ = L '°fog3 + 1 and split the sum above into terms with j' < j and the 
rest with j' > j, where j = d + n — 2 + £. We get, 

fc-2 2j' _ ~i °° „2j' 

V - < V T{d+n-2-j'){j' + l) , < 

1/(1-2) 

Here, in the first sum, we bounded {d + n — 2 — + 1) (with maximum at j' = d+n-3 ^ -^y 
|((i + n — 1)^ and, in the second sum, we took into account that L > 3 and thus c^L~^ < ^. 

Combining f l4.66p and fl4.69p with (14. 70 p . we get the sought bound for (I4.59P and (I4.6ip with 
constants c' = c(2 + jc^^) and r] = max{^{d + n — lY + 2, d + n + 8) . For (I4.60p . the constants 
must be increased by adding 8 to ?7 and multiplying c' by c. □ 

Proof of Theorem [2711 By (|446|) . Lemma gJl and the bound (1420]) we have ||Cfe(p)|| < 
2cMk-i,cApWk-iMp)\p\~^ + cMk-i,cApWk-iMpM^^ Now, for /c = 1, . . . , iV,^he desired 
bounds follow from Lemma |M] and f09|) . For /c = iV + 1 we use flOTI) to get ||CAr+i(p)|| < 
cMjv,c,i,(p)^ and the assertion follows again from Lemma [4.31 and (14 .49 p . □ 



5. Analytic dependence on A and proof of Theorem 12.21 

We now study the dependence of the finite range decomposition on the map A which appears 
in the operator A = V*A'V. To this end we will show that the operators Ux, T and "Jl can be 
locally extended to holomorphic functions of A and the bounds derived previously for fixed A 
can be extended to a small complex ball. Then the Cauchy integral formula immediately yields 
bounds on all derivatives with respect to A. We do not claim that the extensions of U^, 7, Jl to 
complex A yield a finite range decomposition for complex A (indeed positivity is meaningless 
if A is not Hermitian). The extension is merely a convenient tool to show that the relevant 
quantities are real-analytic as functions of real, symmetric, positive definite A. 

Let A be a linear map from c'"^'' to c™^'^ such that 

A = Ao + Ai (5.1) 

with Aq and Ai satisfying the following conditions 

(^F,G')c-xd = (F,AoG)c-xd, (AoF,F)cmxd > Co|F|^ for all F, G G C"^^ and (5.2) 

ll^ill < |- (5.3) 

Here, Cq > is a fixed constant and, as before, (■, ■)cmxd and |-| denote the standard scalar 
product and norm on c™^'^ and \\Ai\\ is the corresponding operator norm of Ai. 

Again, we consider the operator 

A := V*AV, (5.4) 

on ATjv, i.e., 

d 

{A(p)r := ''^V*{AVip)j^r, where {Vip)j^r = ^jfr- (5.5) 
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With operator A we associate the sesquihnear form 

{ip,^)A:= {AVif.Vij) (5.6) 
where (■, ■) is the ^2-scalar product on X^, defining the adjoint A* by 

(Aif,^) = {ip,ij)A = {v,A*ij), with yi* = VM*V, (5.7) 

where A* is the adjoint of A. Note that for real, symmetric A the form (-, ■)a is a scalar product 
and agrees with (■,■)+• ^^e following, we use the previous notation 71+ for the Hilbert space 

1 /2 

with the scalar product (■, ■)ao and define Hv^lUo •= (v^!V^)ao • 

Using and z* to denote the real part and the complex conjugate of a complex number z, 
we summarize the main properties of the sesquilinear form (-, 

Lemma 5.1. Assume that an operator A satisfies the conditions (15.11) . (15. 2p . and (15.31) . 

Then the sesquilinear form {■,-)a on Xn satisifes 

^{^,^)A>'^y\\X, (5.8) 

|(</',^)a| < illv'lUoll^lUo, (5.9) 
{i;,<f)A = {f,i^rA*. (5.10) 

Proof. The first claim follows using the definition of the form (-, ■)a and the lower bound 

^{AF, F)cmxd > {AoF, F)c™xd - j\Ff > ^{AqF, F)c™xd (5.11) 

implied by (Q and f lOl) . 

Using (15. 3p . the Cauchy-Schwarz inequality for the scalar product {AoF,G)Qmxd, and the 
bound from (15. 2p . we also get 

\{AF,G)cmxd\ < {AoF,G)c^xd + -^\F\\G\ < {AqF, F)l''^^a{AQG,G)l^2xd+ 

2 (5.12) 

implying the second claim. 

The last identity follows from the relation {AG, F)cmxd = {G, A*F)cr^xd = {A*F, G')*„xd- □ 

In view of the above Lemma, the complex version of the Lax-Milgram theorem can be used 
to ensure the existence of of the bounded inverse operator = A~^. 

In the following, similarly as in the case of the Hilbert space we use Hj^{Q + x) to 
denote the corresponding Hilbert space (of functions from Xjsi with support inQ + x) with the 
scalar product (■, Oaq- 

Next, we define an extension of the operators II^ for a general complex A. 

Lemma 5.2. Assume that A satisfies (15. ip . (15. 2p . and (15.30 . Then, for each ip ^ X^q, there 
exists a unique v G 1-L+{Q + x) such that 

{v,i,)A = {^,i^)A for all ^en+iQ + x). (5.13) 



Proof. The assertion follows from Lemma 15.11 and the Lax-Milgram theorem. 



□ 
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Lemma 5.3. Assume that A satisfies (15. ip . (15.21) . and (15. 3p . For any ip G Xn, we set 

nA,x^P ■= V, Ha ■■= IlAfi, (5.14) 

with V G 1-L+{Q + x) defined by (I5.13p . Using, as before, to denote the translation by x, 1q 
for the characteristic function of a set Q, and D for the open unit disc D = {w & C : \w\ < 1} , 
we have 

(i) IIa^xTx^ = T^nAf, 

(ii) iiiiAV'lUo < ^f\\Ao, 

(iii) IlAf = (p for all ip G H+iQ), UaUa = Ua, 

(iv) TJaV? = (filq if<f = on Q\Q, 

(v) {Ha^P, ij)A = (</3, nA-ij)A, 

(vi) The map z H- iT^^+^^^y? is holomorphic for z in the open unit disc D. 
Proof. 

(i) : Given that the shift T-x is an isometry with respect to {■■,-)ao ^"^^ maps 'H+{Q + x) onto 
l-L+iQ), we have the identities {TxnA'P,ip)A = {nA'P,T_xi^)A = {'P,T-x4^)a = {Tx<p,i')A for all 
Ip G l-L+iQ + x). As t^Ua'P G H+i^Q + x), this yields the assertion. 

(ii) : Taking ip = varPiA'P in the definition (15.130 of f = IIa'P and using Lemma [5.11 we get 

WWa^WX < ^{nA^,nA^)A = ^{p>,nA^)A < iMUAnA^W (5.i5) 

(iii) : The second assertion follows from the first. By definition, we have {IIa'P ~ fi'ip)A = 
for all Ip G 1-Lj^{Q). In particular, by the assumption ip G 'H+{Q) we may take ip = IIa'P — f 
inferring that IIa'P = 'P- 

(iv) : Let ip> = p>1q. By (iii) we have IIa'P = 'P- Moreover p> — 'p vanishes in Q. Thus V{'p — 'p) 
vanishes in Q-. Hence {'p — 'P,iP)a = for all ip G l-L+^Q) since Vip is supported in Q_. 
Therefore IIa{'P — <^) = which yields the assertion. 

(v) Since IIa'^J G 1-L^{Q) we have 

((/?, nA'^)A = (iIa^, iIa'-iP)a = {iIa*iP, iIa^Ya* = (V^' nA'pYA, = (nA^, ^p)A, (5.16) 

where we used the relation {(p,ip)A = {ip,(p)% and the definition of II a* ■ 

(vi) : This follows from the complex inverse function theorem. Fix ip and consider the map R 
from D X H^{Q) into the dual of 1-L^{Q) given by 

R{z,v){iP) = {v-ip,iP)Ao+zA, = {iAo + zAi){\/v-\/ip),\/iP). (5.17) 

Then R is complex linear in z and v and hence complex differentiable. By the definition of IIa 
we have R{z, w) = if and only iiv = IIaq+zA^'P- Finally the derivative of R with respect to the 
second argument is given by the map Iz from l-L^i^Q) into its dual with Iz{v){ip) = {v, iP)ao+zAi- 
By the Lax-Milgram theorem, 1^ is invertible for z E D. Hence the map z i— ?■ IlAg+zAif is 
complex differentiable in 2;. □ 

Note that for a real symmetric A the above definition of IIa,x agrees with the definition of 
Tlx in Section [3l We define, as before, 

7a ■■= Yl = id - 7a- (5.18) 

Then the following weaker version of Lemma 13.41 holds. 
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Lemma 5.4. Assume that A satisfies (15. ip . (15.21) . and (15. 3p . Then 

\\7a^\\ao < 9||vlUo for all ^ e Xn- (5.19) 

Proof. This is an adaptation of the argument from |1] to the complex case. For the conve- 
nience, we include the details. We have 

I^'W^a^Wao <2f^\i7A^,7Aip)A\ < 2 Yl \inA,.^,nA,y^)A\. (5.20) 

x,j/eTjv 

Set := V TIa^xV- Then vanishes outside Q_ +x since TIa^xV vanishes outside Q + x. Thus, 
in view of (15. 6p and (15. 9p . we get, similarly as in (I5.12p . 

\{nA,x'^,nA,y^)A\ = \ {AT^,Ty) \ = \{AlQ^+^T.j.,lQ_+yTy)\ = \ {AIq _ +yT^ , 1q _ +^Ty) \ < 

< |(y4ollQ_+j;T2:, llQ_+yTj.)"^^^(Ao]lQ_+a:T'y, Ig^+a-Ty)"'^/^ < 

< z{^o'^Q-+yTx, "^Q^+yTx) + | (Aolg.+xT'y, Hq.+i^j/) = 
= l{^O^Q-+yTx,Tx) + ^{AolQ__^.^Ty,Ty). 

(5.21) 

Now Xlyeijv ^Q+y constant function l"^ and thus 

5^ I {Ha^x^, nA,yip)A\ < Yl ^"^oTx, Tx) = li" inA,x^, nA,x^)Ao < 

= 3l'''^i^,7Aif)A < II'^MaoW^a^Wao- 

Combined with (I5.20p . this yields the assertion. □ 

Next, we bound the Fourier multipliers of operators 7a and JIa- Using the relation Ua^x = 
TxUaT-x we get, as before, 

7Aiafp) = I-'' Y e-"^'"'^nA{aU){z) U (5.22) 

and thus, the Fourier multiplier Ta{p) is given by 

fA{v)a = l-" Y e"*^^''^^A(a/p)(z). (5.23) 

zGTjv 

Also, the operator A = V*AV satisfies again the equation 

A{afp) = {A{p)a)fp (5.24) 

with 

{A{p)a, b)cm = {A{a ® q{p)), b ® q{p))c^xd, q{p)j = e'^^ - 1. (5.25) 

Hence, 

\\Ap)\\ < fPoll \p\' and WAipr'W < (5.26) 

Lemma 5.5. Assume that A satisfies (15. ip . ( 15. 2p . and (15. 3p . Then, there is a constant c < oo 
(depending only on ||ylo||,Co, and d) such that, for all p G T^r \ {0}, 

(i) \\Mp)'^'fA{p)Mpr'^'\\ < 9, \\Mpy^'T^Aip)Mpr'^'\\ < lo. 



FINITE RANGE DECOMPOSITION 27 



(ii) \\fAip)\\<cmm{l,i\p\l)^). 

(iii) <cmin(l,l(l + ^)). 

Proof. The proof is largely parallel to the proof of Lemma 1411] for real symmetric A, but the 
constants are slightly worse since {■,-)a is no longer a scalar product. 

(i) : The second bound follows from the first since 1Za{p) = ^ — Ta{p)- For the first estimate, 
we apply Lemma with (p = afp, a G C™. This yields 

{Ao{p)fA{p)a,fA{p)a)c^^ < 81{Ao{p)a,a)c^. (5.27) 

Taking a = Ao{p)^^^^b, we deduce that 

(A(p)'/'fA(p)lo(p)-'/\ A(p)^/'f^(p)lo(p)-'/'&)c™ < 81{b,b)c^ (5.28) 
and this finishes the proof of (i). 

(ii) : First, we assume that \p\l > 1. It follows from (i) that 117^(^)11 < 9|| A(p)~^/^|| II A(p)^^^||- 
By g;20D we have ||A(p)'/'|| < Po||'/'bl and WMpy'^'W < ^^^p, yielding (ii). 

Now, assume that \p\l < 1. We first estimate the norm of -4.(p)7a(p)- Note that 

l''{Aip)fA{p)a,b)c^ = {A{p) J2 nA{afp){z),bfp{z))c^ = {nA{afp),A*{bfp)) = {nA{aU)Mp)A. 

z&Tn 

(5.29) 

Thus by Lemma [5l3] (iii), (v), and (15.91) . we get 

l''\{A{p)fA{p)a,b)c^\ = \{nAnA{afp),bfp)A\ = \{nA{afp),nA4bfp))A\ 

<f||i7^(a/,)|UJ|i7^.(6/,)|U,. 



(5.30) 



To estimate ||i7^(a/p)||Ao we use (15.81) and the fact that (IlAi^p), nA{^p))A = {'^ , n a{'^)) a 
according to (I5.13P since TIa^^) G H+{Q), yielding 

|||i7^(a/,)||^„ < \{nA{afp),nA{afp))A\ = \{afp, nA{afp))A\ = | (yi(a/p), i7^(a/,)) | 

= \{iAip)a)fp,nA{afp))\ = \{iAip)a)ifp-l),nAiafp))\. (5.31) 

In the last step we used the fact that functions in have average zero. Now i7A(a/p) is 
supported in Q and \ fp{z) — 1| < v^|p|/ for z ^ Q. In combination with fl5.26p this yields 

UnAiamX < ||(l(p)a)(/,-l)lQ|| ||i7^(a/,)|| < |Po|| Ip^Td |a| (|p|0/'^/^||i7A(a/,) ||. (5.32) 
The Poincare inequality ( [9l [10] ) implies that 

\\nAiafp)y^Q) < cl\\VnAiafp)y^Q) < ^/||i7^(a/,)|U„. (5.33) 
Combining the inequalities above we get 

||i7A(a/p)|Uo < 3Vdc^^\p\H'l'/'\a\. (5.34) 



Co 



The same estimate holds for Ua*- Hence, from ( ]5.30p . 



l^\{A{p)fA{p)a,b)c^\ < fd^^^^lpflWlaWb]. (5.35) 



Co 
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With the help of f l5.26p . this yields the claim for a suitable constant c. 

(iii): For \p\l < 1 the estimate follows from (i). Thus, we assume \p\l > 1. Again we first 
estimate A{p)71a{p) ■ Since 7^a(p) = 1 — Ta{p) we get from fl5.29p 

l''{Aip)nAip)a,b)cm = l''{Aip)a,b)c^ - {nA{aU)M'v)A. (5.36) 
Let be a cut-off function such that 
uj{z) = lii z eQ\Q, w = Oif dist(z,Q\g) > 1 + ^, < u; < 1, |Va;| < c|p|. (5.37) 

By Lemma [573] (iii), we have -//^(l — ^){o,fp) = (1 ~ ^)^Q(^fp- Hence 

{nAiafp),bfp)A = (nAiaufp), bfp)A + (a(l - w)1q/p, bfp)A (5.38) 

and 

ia{l-u)lQfp,bfp)A = {ail-u)lQfp,A*{bfp))A={ail-u)lQfp,iAiprb)fp) 

= ^(l-c^)(l(p)a,6)c".= 5^(l-c^)(l(p)a,&)c".. (5.39) 

In the last step we used that w = 1 on Q- \ Q- Since \Q- \ = l"^, this yields 

l''{A{p)1ZA{p)a, b)cm = -{nAiaoofp), bfp)A + uj{A{p)a, b)^^. (5.40) 

Given that u is supported in a neighbourhood of the order 1 + l/|p| around Q \ Q, the last 
term is bounded by 

I 1^ u{A{p)a,b),..\<4dl'^\l + ^^)\\Aip)M^^^ (5.41) 

To estimate the remaining term we introduce again an additional cut-off function a) for which 

2 3 

cl;(^) = 1 if dist(2,Q\Q) < 2+ — , w = if dist(^, Q\Q) > 3+ — , < w < 1, < c\p\. 

\p\ \p\ 

(5.42) 

Then, taking into account that 1 — u and A{aujfp) have disjoint support, 

{nA{aufp),bfp)A = (aufp, nA*{bfp))A = (aufp, nA-{bujfp))A + (aw/p, (1 - uj)lQbfp)A 
= {aufp, nA*{bujfp))A + {A{aujfp), (1 - co)lQbfp) 

= {aujfp,nA*{bCjfp))A. (5.43) 

Hence, 

\{nA{aojfp),bfp)A\ < f ||acj/p|Uo ||i7A- (&f^/p) IUo < ^aufp\\Ao II^^/pIUo 

< C\\Ao\\\p\Y-\l + ^)\a\\b\, (5.44) 

where we used that u and u are supported in a strip of order 1 -|- l/\p\ around Q \Q, that 
l/bl < ^ and that the gradients of w, and fp are bounded by c\p\. The combination of (15.40 p . 
(I5.4ip and fl5.44p now yields the estimate 

\\A{p)nA{p)\\<C\\Ao\\\prj{l + ^^). (5.45) 

In view of f l5.26p this finishes the proof of (iii). □ 
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We now define and estimate the operators GA,k- Assuming again that L > 16, we use 

Q, = {1, ...,/, - l}'^ with I, = +1 for J = 1, . . . , iV, (5.46) 

to define Taj,T^j-, "J^aj, and ^'aj, in the same way as before. Introducing 

MA,k:=^A,i...^A,k-i^A,k, for k = l,...,N, and Ma,o := id, (5.47) 

we set 

eA,k := MA,fc-ie^M'^_,_, - MA,kGAM'^^i,, (5.48) 

for /c = 1 , . . . , and and 

6^,^+1 := M^.^CaM'^^^ (5.49) 

for = AT + 1. 

Considering the same annuh Aj, j = 0,1, N, introduced in fl4.50p and f l4.5ip . as well as 
the functions Mk^chip) and Mk^ciip) defined in fl4.52p and fl4.53p . we have the following bound. 

Lemma 5.6. Assume that A satisfies (15.11) . (15. 2p . and (15. 3p . Then there exists a constant 

c (depending only on ||y4o||,Co and the dimension d) such that for all L > 16, all N, and all 
j = 1, . . . N , we have 

\\MaAp)\\ <cMk,cAp) (5-50) 

for k = 0, . . . , N and 

\\fA,k+i{p)MAAp)\\ <cMu,cAp) (5-51) 

fork = 0,...,N~l. 

Proof. This is similar to Lemma [4.21 for the case of real symmetric A. However, the bound 
||'7^(p)|| < 1 for < j has to be replaced by 

\\nA,k{p)\\ < 1 + \\fAAp)\\ < 1 + c{\p\L'')\ (5.52) 

This yields 

j j 

n W^aAp^ < + c{\p\L'Y) < c, (5.53) 

k=0 k=0 

resulting in the additional factor c in the claim. □ 

With the help of Lemma [4.31 we can now bound ||Cyi,fe|| in the same way as in the proof of 
Theorem 12.11 starting from (I5.49P for = A^ + 1 and from the equality 

Cyl.fc = (X4,feMyi_fc_i)CAM^ + MA,fc-lCA('^A,fcMA,fc-l)' — 'J'A,A:MA_fc_ (5.54) 

for k = 1, . . . , N. The latter follows from (I5.48P using the equation 3^a,A; = id — 7A,k- 

We can now use the Cauchy integral formula to control the derivatives with respect to A. 

Lemma 5.7. Let D = {z e C : \z\ < 1}. 

(i) Suppose that f : D ^ Qmxm holomorphic and 

sup 11/(2)11 <M. (5.55) 

Then the j-th derivative satisfies 

11/(^^0)11 <Mj! . (5.56) 
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(ii) Suppose that f : D ^ Qmxm g - d ^ ^rnxm holomorphic and 

sup 11/(^)11 < Ml, sup \\g{z)\\ < M,. (5.57) 
zeD zeD 

Then the function h(t) = f{t)g*{t) is real- analytic in (—1, 1) and 

||/i(-')(0)|| < MiMaj! . (5.58) 
Here g*{t) denotes the adjoint matrix of g{t). 

Proof. Assertion (i) follows directly from the Cauchy integral formula. To show (ii), we note 

that g{z) = ^jajz^ with aj G C"^^". Define G{z) := EjS*^^- ^hen G{z) = g{z*)*. Hence 

||G(2;)|| = ||5f(2;*)||. Thus H := fG is holomorphic in D and satisfies sup^, \\H\\ < M1M2. Hence 

i7(^)(0) < k\ M1M2. For t E (-1, 1) we have H{t) = h{t) and the assertion follows. □ 

Proof of Theorem 12.21 It only remains to show the claim (iii). Let Aq and Ai as before 
and assume in addition that Aq and Ai are real and symmetric. Set 

A{z) := Aq + zAi (5.59) 

Then the maps 

Z ^ Ma{z)Ap)^ Z ^fA{z),k{p)MA(z),k{p), Z ^ MA(z),k{p)CA{z){p), 

Z ^ T'A{z),k{p)MA(z)M^A{z){p) (5.60) 

are holomorphic in D. Moreover 

\\Ma(z)Ap)\\ < cMk,cAp)^ 

\\fAiz)AP)^Aiz)AP)\\ < cil4,c,L(p), 



\\MA{z),k{p)^A{z)ip)\\ < j^Mk,cAP), 



\\fA(z)Ap)MAiz),k{p)CA(z){p)\\ < j^Mk,cAp), (5.61) 

again with the functions M^^cl and M^^cl defined in fl4.52p and fl4.53p . respectively. Hence it 
follows from Lemma [5.71 that 



d^' 



dp \t=o 



cjl 



— CA,+tA,Ap) <T^{'2Mk-i,cAp)Mk-i,cAp) + M^-i,c,Ap)) (5.62) 



\P\ 



for k = 1, . . . ,N (with the obvious modification for k = N + 1). Thus Lemma 14.31 and the 
estimate fl4.49p for the inverse Fourier transform yield 

sup \\{V''D{CaoAx){Au . . . ,Ai)\\ < C„(d)jU-('=-i)('^-2+l°l)L''("''^). (5.63) 
Finally suppose that \\A\\ < 1 and set Aq = ^A. Then the desired estimate f l2.24p follows from 



FINITE RANGE DECOMPOSITION 



31 



References 

[1] R.A. Adams and J.J.F. Fournier, Soholev Spaces, Academic Press 2nd ed., Elsevier, (2005) 

[2] S. Adams, R. Kotecky and S. Muller, Strict convexity of the surface tension for non-convex 

potentials, in preparation. 
[3] D.C. Brydges, Lectures on Renormalisation group, in Statistical Mechanics, editors S. Sheffield 

and T. Spencer, IAS/Park City Mathematics Series Vol. 16, American Mathematical Society, pages 

7-93, (2009). 

[4] D.C. Brydges and A. Talarczyk, Finite range decompositions of positive-definite functions. 
Journal of Functional Analysis 236, 682-711 (2006). 

[5] D.C. Brydges, G. Guadagni and P.K. Mitter, Finite Range Decomposition of Gaussian Pro- 
cesses, Journal of Statistical Physics 115, Nos. 1/2, 415-449 (2004). 

[6] D.C. Brydges and H.T. Yau, Grad Perturbations of Massless Gaussian Fields, Commun. Math. 

Phys. 129, 351-392 (1990). 
[7] T. FuNAKi and H. Spohn, Motion by Mean Curvature from the Ginzburg-Landau Interface 

Model, Communications in Mathematical Physics 185, 1-36 (1997). 
[8] C. Hainzl and R. Seiringer, General decomposition of radial functions on 'R'^ and applications 

to A-body quantum systems, Lett. Math. Phys. 61, 75-84 (2002). 

[9] B.G. Pachapatte, , On discrete inequalities of the Poincare type, Periodica Mathematica Hun- 
garica 19, 227-233 (1988) 

[10] W.S. Cheung, Some discrete Poincare-type inequalities, International Journal of Mathematics and 

Mathematical Sciences 25, 479-488 (2001). 
[11] K.G. Wilson, The renormalization group and critical phenomena. Rev. Modern Phys. 53 (3), 

583-600, (1983). 



